Magneto-optical conductivity of Weyl semimetals 
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Weyl semimetals are a topological phase of matter that have drawn recent interest and have 
been suggested as a possible phase of the pyrochlore irridates, among other materials. Here, we 
compute the magneto-optical response of an isolated Weyl point. We find that the conductivity is 
a series of asymmetric peaks lying on top of a linear background. We find additional features in 
the conductivity that are explained by the underlying Landau level structure. We show that weak 
disorder tends to blur out the peaks, but they still show up as measurable oscillations on a linear 
background. 
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I. INTRODUCTION 

A Dirac semimetal is a material in which the con- 
duction and valence bands only touch at isolated points 
within the Brillouin Zone (BZ). Around these points, the 
dispersion is linear and the low energy theory is described 
by a Dirac Hamiltonian. Perhaps the most famous Dirac 
material is graphene, which contains 2 inequivalent Dirac 
points in its BZ. These two dimensional Dirac points are 
not robust against perturbations. Any perturbation pro- 
portional to fJz will gap out the band touching. For ex- 
ample, in graphene spin orbit coupling should induce a 
mass gap. However, since carbon has such a small atomic 
number, spin orbit coupling is very weak and this split- 
ting has never been measured in graphene. 

Dirac semimetals can also exist in three dimensions. 
An isolated band touching is described by the Hamilto- 
nian 



T-i — vpcr ■ k. 



(1) 



where vp is the Fermi velocity, cr is the vector of Pauli 
matrices, and k is the momentum as measured from the 
band touching. This equation is the electronic analogue 
of the Weyl equation from particle physics. Materials 
described by such a Hamiltonian have become known as 
Weyl semimetals. The Weyl Hamiltonian is completely 
robust against any perturbation since it uses all three of 
the Pauli matrices. In materials with time reversal sym- 
metry (T) and inversion symmetry (I) then there must 
be four bands linearly dispersing around any band touch- 
ing point in the BZ. Such a material is called a 3D Dirac 
semimetal and it is not robust against perturbations since 
there are additional Dirac matrices in the 4x4 represen- 
tation. In a system with broken T or I it is possible to 
have a phase described by the Weyl Hamiltonian. 

Weyl points also have interesting topological proper- 
ties. They are monopoles of Berry flux in the BZ. This 
peculiar momentum space topology manifests itself in 



several ways. One way is the peculiar nature of the 
gapless surface states. The surface states of a Weyl 
semimetal are known as Fermi arcs.^ These Fermi arcs 
are open segments of Fermi surface connecting the projec- 
tions of the Weyl points onto the 2D surface BZ. Another 
manifestation of the topological nature of Weyl points is 
in the low energy electromagnetic response. The low en- 
ergy effective action contains a topological 6'-term and is 
described by axion electrodynamics.^ The physical effect 
of this exotic 0-term is an anomalous Hall conductivity;^ 
The equations of motion also show that these materi- 
als could exhibit a chiral magnetic effect ~ dissipation- 
less currents in response to an applied magnetic field. 
Recent numerical experiments found no evidence for the 
chiral magnetic effect using their lattice model for a Weyl 
semimetal.^ Charge transport in Weyl semimetals in the 
presence of disorder have also been studied, where they 
found two distinct regimes for the conductivityi^ The di- 
rect observation of Fermi arcs or the anomalous Hall ef- 
fect would be good evidence for the existence of a Weyl 
semimetal. 

Although Weyl semimetals have not yet been ob- 
served there are several candidates. The pyrochlore 
irridateai*^^— are a promising host for the Weyl semimet- 
alic phase. Proposals have also been made for the ex- 
istence of a Weyl semimetalic phase in topological in- 
sulator heterostructureSf^"— as well as in magnetically 
doped topological insulators. ^^■^'^ A recent paper shows 
that it may even be possible to induce a topological phase 
transition to a Weyl semimetal using the electron-phonon 
interaction.-" 

The optical conductivity in a variety of quasicryst aline 
materials is linear over a broad range in frequencies!^ 
This has been interpreted as evidence for Dirac fermions 
in these materials. Further evidence is needed to show 
that the physics in the quasicrystaline samples is due to 
Weyl fermions and not simply Dirac fermions. 

In this paper we compute the magneto-optical conduc- 
tivity for a single Weyl point. For N Weyl nodes the 
conductivity simply scales by a factor of N. We find that 



the transverse optical conductivity is a series of asym- 
metric peaks with onsets proportional to the square root 
of the magnetic field. These peaks lie on the background 
of the zero magnetic field limit. The structure of these 
peaks as the chemical potential is changed is rich, and its 
behaviour follows from the underlying dispersive Landau 
levels. We obtain formulae for the absorptive parts of 
both the transverse and Hall conductivities in the clean 
limit. Lastly, we examine the conductivity for circularly 
polarized light, as well as the semi-classical limit. These 
all provide verifiable predictions distinct to Dirac physics 
that can give additional support for the existence of a 
Weyl semimetal. 



II. DYNAMICAL CONDUCTIVITY OF A WEYL 
SEMIMETAL 

The Hamiltonian for an isolated Weyl point in a mag- 
netic field is given by 



H 



(T, 



(2) 



where e is the electron charge, and c is the speed of light. 
In the above we have ignored the effect of Zeeman split- 
ting, since this effect should be small at accessible mag- 
netic field strengths. We consider a field B = Bz, and 
use the gauge Ay = — Q and A^ = —By. We work in 
units where h = ks = vp = ^- It is convenient to define 
the magnetic length which is given by = c/eB. We 
define the operators 



V2 



(3) 
(4) 



where 11 = — iV H- eA. Then, the Hamiltonian is given 
by the following set of 2 x 2 matrices 



n = 




(5) 




FIG. 1: Dispersive Landau level structure for a Weyl 
semimetal at finite chemical potential. In this figure the chem- 
ical potential lies between the fifth and sixth Landau levels. 
The arrow at fcz = shows an interband transition from n = 5 
to n = 6. These interband transitions at fcz = give the main 
set of peaks shown in Figure [21 The two transitions at finite 
occur between the n = and n = 1 Landau levels. They 
occur when ^ lies between two Landau levels and are respon- 
sible for the additional bumps seen in Figure [4] 



This is the same equation as for a superconductor, with 
the identification = e, En\ = E and \/2n/lB = A. So 



VnX = 




(8) 
(9) 



First we find solutions to the Eigenvalue problem H^p 

VnX 



E^p. We take ?A = 



and obtain 



En = ±4 



I2n 



and 



VnX- 



k'i — Xt 



I2n 



2n 



k-' = E, 



71 X ■ 



UnxMEnX - fcz). 



(6) 



(7) 



These equations describe the Energy level structure 
and wavefunctions for the set of 3D Landau levels gener- 
ated by the magnetic field. Unlike in the 2D case these 
Landau levels are dispersive (in kz) as shown in Figure 
[T] The structure of these energy levels will control the 
shape of the dynamical conductivity, as we will elaborate 
on later. 

The dynamical conductivity tensor can be obtained 
from the Kubo formula. Expressed in the Landau level 
basis in the clean limit we have 
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FIG. 2: (Color online) The real part of the optical conductiv- 
ity (in units of /SttIb) for fils ~ 0,1.4, and 2 as a function 
of ijjIb- In this figure we took TIb ~ 0.01. The conductivity 
is a series of peaks which sit on a linear background. These 
peaks occur at the kz = value of the Landau levels shown in 
Figure[T] The inset shows the n = conductivity over a larger 
frequency range, demonstrating its linearity. The background 
conductivity should be the result for free fermions, which is a 
straight line of slope 1/3 in our units. 
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2-K V E„ 



Eni) 



OJ + EnX- En'X' + iO+ 



(10) 



where f{x) = 1/(1 + e'^'^"-^)), (3 is the inverse temperature, and /i is the chemical potential. For the dissipative 
components of the conductivity tensor, '^{(t^x) and '^{<Jxy) we obtain 



( f{E^) - f{En+l)\ ({En - k,){En+i + k,) 



V En — En+l J V EnEn+1 

f{En+l) - fi-En) + f{En) - f{-En+l) \ P 



En + En+l J V EnEn+1 

fi-En+l) - fi-En)\ {{En+l - fc.)(S„ + fc,) 



6{UJ + En - En+l) 

5{lO — En — En+l] 



En — En+l 



EnE. 



n-t^n+1 



5{UJ - En+l + En) 



(11) 



^ / 2tt 



dkz 



ff{En+l) - f{En)\ ({En - k,){En+l + fc.) 



V En — En+l J V EnEn+1 

f{~En)-f{En+l) + f{En)-f{-En+l)\f. ^ 



En + En+l 



S{UJ + En - En+l) 

S{U - En - En+l) 



EnE. 



'n+1 



.f{-En+l) - f{-En)\ ({En+l ^ k,){En + fc,) 



En — En+l 



EnEn^ 



6{UJ - En+l + En) 



(12) 



r 



In this problem the magnetic length, I b sets a natural This choice of units will allow us to produce a single plot 
energy scale for the problem. So, for convenience, we for all magnetic fields. In the clean limit we can use the 
define u; = Iguj, k = Igk^, and £■„ = l^En = y/2n + k'^. delta functions to do the integration over k^. We find 
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that there are contributions from and we obtain the following formulae for the absorptive 

part of the conductivity tensor: 

- - 1 /4-4(:;2(2n + l)+w4 



sinh( 



2+1^ 
20T 



COSh(f;)+COSh(^) 



sinh(- 



1 i^^J 



2aT 



cosh(f ) +cosh(^) 



E 

ri=0 



|2(2n+ 1) 



ai^cj4-4cj2(2n + l) + 4 



(14) 



e'^/'^ + cosh(' 



2S}T 



+ cosh(^) 



cosh(f ) +cosh(^ 



cosh(f)+cosh(^) 



E 



|2(2n+l) 



(|V2- 



- w 



(15) 



- nL = 0.0 - 




FIG. 3: (Color online) The real part of the optical conduc- 
tivity (in units of e'^/STr/s) for graphene in the clean limit. 
Here TIb = 0.01. For = there is a series of peaks at 
cj = •\/2(n + 1) + y/2n. As /i is increased these peaks disap- 
pear and are transfered into the intraband peak at low fre- 
quency. 



In the above expressions [xj denotes the integer part 
of X. '^{(Txx) contains a series of peaks that sit on a linear 
background (Figure [5]) . The background is the result for 
free Dirac fermions in absence of a magnetic field, namely 
(To = ttiG units of our plots this is a straight 

line with slope 1/3. In the inset of Figure we show the 
real part of Uxx over a larger frequency range, showing 
that it becomes linear at large Q with small oscillations 
from the Lan dau level structure. For = the peaks oc- 
cur at cD = •\/2(n -I- 1) + \/2n for integer n. The series of 



peaks corresponds to allowed interband transitions in the 
Landau level structure. This peak spacing is proportional 
to \/B. From an experimental point of view this means 
one can see well spaced peaks even for modest fields. The 
asymmetry of the peaks is refiected by the square root 
singularity in Eq. 1141 physically it comes from the disper- 
sive structure of the Landau levels. The long tails of the 
peaks originate from the square root singularity and add 
together to provide the linear background we observe. In 
Figure [5] we have shown the real conductivity for 3 values 
of the chemical potential. These values of the chemical 
potential are given by the energies of the Landau levels. 
We see that as the chemical potential is increased, the 
peaks at low energies disappear and the optical spectral 
weight is entirely transfered into an intraband peak at 
low frequency. The pattern for how the peaks disappear 
is reminiscent of the what happens in grapheneJ^ For 
comparison we have calculated the optical conductivity 
for graphene (Figure [3|). Indeed, the peaks occur at the 
same positions as for a Weyl semimetal, and disappear as 
the chemical potential is increased. However, in graphene 
the Landau levels are non-dispersive, and changes to the 
optical conductivity can only happen when the chemi- 
cal potential passes through a Landau level. In a Weyl 
semimetal, the dispersive nature of the Landau levels pro- 
duces a richer structure. 

Figure |4] shows how the dispersive Landau level struc- 
ture affects the optical conductivity. The top panel shows 
the evolution as the chemical potential is moved between 
two adjacent Landau levels. Unlike in graphene, the 
spectral weight of the interband peaks is continuously re- 
distributed into the intraband peak. The bottom panel 
shows how the conductivity changes as the chemical po- 
tential moves through a Landau level. As the chemical 
potential passes through a Landau level there is a large 
redistribution of spectral weight, much like in graphene. 
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FIG. 4: (Color online) The real part of the conductivity (in 
units of /SttIb) as a function of frequency for several val- 
ues of the chemical potential. The top panel illustrates how 
the conductivity changes as the chemical potential sweeps be- 
tween adjacent Landau levels. The bottom panel illustrates 
the change in the optical conductivity as a Landau level is 
crossed. 



The lowest peak in the sohd black curve disappears en- 
tirely as ^ sweeps through the energy of the first Landau 
level, while the rest of the curve is unaffected. It reap- 
pears as a peak in the dashed orange curve that comes 
from the intraband transition between the n = 1 and 
n = 2 Landau levels. The tails of this new peak overlap 
with the original curve. As fi is incremented further to 2 
(the energy of the n — 2 Landau level) the second peak 
in the black curve also becomes modified and now car- 
ries only half its original spectral weight. In a next step 
the remainder of this second peak will entirely diappear 
into the intraband transitions. The dispersive Landau 
levels decorate the graph with additional features. Both 
the top and bottom panel of Figure H] show additional 
shoulders that are not located at the signature frequen- 




FIG. 5: The real part of the optical conductivity (in units of 
e'^ /8tvIb) for left-handed circularly polarized light for fl — 2.0 
and T = 0.01. In the low temperature limit, the imaginary 
part of the Hall conductivity exactly cancels the real part of 
the diagonal conductivity for uj < 2/i. 



cies w = •\/2(n + 1) + y/2n. These additional features 
are interband transitions involving the zeroth Landau 
level. More precisely they are the n — n = 1 and 
n = 1 — > n = transitions. Since the zeroth Landau 
level is an envelope for all the other Landau levels these 
additional features exist for all values of the chemical po- 
tential. For a chemical potential fl the shoulders onset at 
u) = fi+ and ui = p, + \JjiF—2. 

So far we have only discussed the diagonal conductiv- 
ity "^(axx)- For experiments that probe the polarization 
of light, such as the Faraday and Kerr effects, the quan- 
tity G± = Gxx ± ic^xy is appropriate. (T+ describes light 
with right handed polarization, and (T_ describes light 
with left handed polarization. The absorptive part of g± 
is 3fi((T±) — '^{oxx) T S(crj,y). A plot of 3?(o-_) is shown 
in Figure [5] For left handed light there is a cancella- 
tion between the Hall conductivity and the longitudinal 
conductivity in the frequency range w < 2[i. For right 
handed polarization (not shown) the peaks in the region 
a; < 2/i have double the magnitude as in Figure [2] The 
presence of peaks for left circularly polarized light in this 
region could be a signature of correlation effects. 

To realistically treat the problem of disorder (espe- 
cially at large scattering rates) one would need the study 
of the Kubo formula expressed in terms of the spectral 
density, A(k,iS). By solving a Boltzmann equation, one 
could obtain an expression for the self energy. The imag- 
inary part of the self energy then serves to broaden the 
spectral densities. Here we capture the effect phenomeno- 
logically by introducing a scattering rate, F. By making 
the replacement 



<5(z)^ 



1 



TT r2 + z2 ' 



(16) 
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FIG. 6: (Color online) The real part of the optical conductiv- 
ity (in units of e'^ /SttIb) for ft = 0.0 and T — 0.01 for several 
values of scattering rate, f . The red line shows the conduc- 
tivity in the absence of a magnetic field. Increasing disorder 
tends to blur out the peaks towards the free limit. 
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FIG. 7; (Color online) The real part of the optical conduc- 
tivity (in units of /S-kIb) in the semi-classical limit. This 
occurs for large values of the chemical potential. All the spec- 
tral weight of the now forbidden interband transitions piles up 
into a single resonance at low frequency. 



in Eqn [TT] we can understand how weak scattering will 
effect the conductivity. This replacement should give a 
qualitatively correct picture for small values of the inverse 
scattering time, F. 

Figure [H] shows the real part of the transverse con- 
ductivity for f = 0.05, 0.1, and 0.2. The straight red 
line shows the conductivity of free 3D Fermions, cto for 
comparison. The disorder tends to smear out the peaks, 
and at large wl b the conductivity trends towards the free 
limit as the disorder is increased. In particular the scat- 



tering rate needs to be large compared to the energy level 
spacing to see any peaks. In Figured the first few peaks 
are still very well defined for f — 0.2. For a Fermi ve- 
locity of vp ^ lO^m/s, and B measured in Tesla, this 
corresponds to a scattering rate of F '-^ 7\fBmeV . 

Finally we consider the semi-classical limit. The semi- 
classical limit occurs when the Landau level quantization 
is no longer important. This occurs when the chemical 
potential is the largest energy scale in the problem. For 
i?jv < p. < En+1, with ^ 1 we have that 

= En+i - ^at ~ 4. (17) 
M 

In this limit, if one goes back to the unbarred units, the 
energy spacing goes like B in contrast to the ^/B be- 
haviour at small ^. The semi-classical resonance it shown 
in Figure [7] for /2 = 5, and 10. It appears as a peak at 
l//i, and contains the spectral weight from the interband 
transitions that have disappeared between Q < uj < 2/1. 
In graphene where the Landau levels are flat this line is 
due only to the Ej^ — >■ E^^i intraband transition. In the 
present case however, the dispersive Landau level struc- 
ture gives contributions to this lineshape from the set of 
transitions i?„ — >■ En+i with n € [0, A^]. We can also see 
in Figure [7] that as we doubled /x the amount of area in 
the semi-classical peak quadrupled. This is a reflection of 
the linear background of the conductivity. For a chemical 
potential of /i, interband transitions are forbidden out to 
Lu = 2/i. The linearity of the conductivity combined with 
a sum rule shows that the area under this peak should 
scale like jj?, as we observe. 

III. DISCUSSION AND CONCLUSIONS 

In this paper we calculated the optical conductivity 
for a Weyl semimetal in the presence of a magnetic field. 
We found that the dispersive Landau level structure pro- 
duced a set of asymmetric peaks on a linear background. 
The peaks have a ^fB spacing, a signature of Dirac 
physics. The shape of the peaks reflects a square root sin- 
gularity, and their long tails conspire together to give the 
linear background in this case. The positions of the peaks 
occur at the same frequencies as in graphene, although 
there are additional features when the chemical poten- 
tial falls between two Landau levels. We also showed 
that for left hand polarized light, that there is a cancel- 
lation for frequencies w < 2/Lt, and peaks in this region 
could be associated with interaction effects. We showed 
that the semi-classical resonance has a typical lineshape, 
but that it consists of many intraband transitions and 
has spectral weight thats scale like /j,^. We also showed 
the effect of weak disorder. The sharp peaks tend to blur 
out towards the linear background as the scattering rate 
increased. This may make the peaks hard to see in dirty 
samples. In particular, the scatting rate needs to be less 
than the Landau level spacing for the peaks to still be 
observable. 
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Finally we would like to make a remark on units. 
Throughout the paper we have used the dimensionless 
units where Ib set the energy scale. For convenience we 
include the conversion to 'real' units. For vp measured 
in to/s and B measured in T we have 



(18) 



This conversion will be useful for comparing to experi- 
mental results. 
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